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Abstract. 

The Laser Interferometer Space Antenna (LISA) will open the low-frequency (0.1-100 mHz) part of the gravitational wave 
spectrum to direct observation. Of order 3600 galactic close binary white dwarfs will be individually resolvable in its all-sky 
spectrum, of which a dozen systems are expected to be on the verge of merger, showing the effects of strong tidal heating 
and/or early onset of tidal mass transfer. Optical study of these systems would provide important insights into tidal dissipation 
mechanisms, and internal heating in merging white dwarfs that sets ignition conditions for potential type la supernovae. 
Theoretical modeling and instrumentation programs are needed now to enable a campaign for optical identifications to exploit 
this opportunity. 

Keywords: binary stars, white dwarfs, stellar evolution, supernovae 
PACS: 97.80.-d,97.20.Rp,97.10.Cv,97.60.Bw 

INTRODUCTION 

At this writing, the Laser Interferometer Space Antenna (LISA) is targeted for launch in the mid-2020's. This mission 
will open a new frequency window for gravitational wave detection, one which includes the expected signals of several 
thousand close double white dwarfs in the final stages of their inspiral toward merger. Among these systems, some 
may have sufficient total mass to qualify as Type la supernova progenitors. The purpose of this contribution is to call 
attention to the challenge and importance of identifying optically close double white dwarfs discovered by LISA on 
the verge of merger. 

LISA will consist of three spacecraft 5 x 10 6 km apart, forming a roughly equilateral triangle, in what is basically 
a Michelson interferometer. This formation of spacecraft will orbit the Sun in the ecliptic plane at 1 AU, 20° behind 
the Earth. The three spacecraft themselves define a plane inclined by 60° with respect to the ecliptic plane. That plane 
precesses about the ecliptic pole as the spacecraft orbit the Sun, modulating the antenna pattern as they do so. As 
LISA orbits the Sun, periodic sources (such as double white dwarfs) will therefore be modulated in amplitude by this 
variation in antenna pattern, and in frequency by doppler motion of the spacecraft array with respect to the Sun. These 
two effects combined provide a signature by which sources can be localized. 

LISA will be sensitive to GW frequencies in the range v « 0.1 to 100 mHz. It is designed to achieve a limiting 
gravitational wave dimensionless strain amplitude of order h ps 10~ 23 at signal-to-noise ratio SNR = 5 in one year of 
operation [7, 8]. In the lower- frequency part of this window, galactic close double white dwarfs are expected to provide 
such a plethora of sources that LISA'S sensitivity will effectively be confusion-limited by this background. Only very 
nearby sources are likely to be detectable above this background. However, as binary systems inspiral toward merger, 
the accelerating strength of GW losses leads to a characteristic power-law decline in the number of sources per unit 
frequency, N v oc v~ n / 3 . When vN v < c/v@ = 10 4 (v® being the Earth's orbital speed about the Sun), the average 
separation of sources in frequency becomes larger than their characteristic doppler width due to spacecraft orbital 
motion. Individual sources are then resolvable. Webbink and Han [10] estimate that this transition occurs around 
V w 3.6 MHz, and that of order ~ 3600 sources could then be resolved. A more careful treatment of source confusion 
can be found in Bender and Hils [1], yielding a similar estimate of the number of GW sources resolvable in frequency. 



THE GRAVITATIONAL WAVE SIGNATURE OF CLOSE BINARY WHITE DWARFS 

The keys to identifying close double white dwarfs on the verge of merger are contained within their gravitational wave 
signatures. Systems at this evolutionary stage are expected to have circular orbits, as universally observed among 
known close double white dwarf binaries. The gravitational wave signal then immediately fixes the orbital period and 
phase of the source (modulo a 180° ambiguity in phase because the GW frequency is twice the orbital frequency: 



^orb = 4n/v). The dimensionless strain amplitude of a source a distance d from the observer is [cf. 6, 4] 
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where m is the "chirp" mass, a function of the binary component masses M\ and M2, 
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and h + and h x and the separate dimensionless amplitudes of the two polarizations of the (quadrupole) gravitational 
wave. The ratio of h x to h+ is a function of the orbital inclination ;' of the binary's orbital axis to the line of sight to 
the observer: 

h x 2 cos; 
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Thus, the gravitational wave signature of a source specifies not only the orbital ephemeris of that source, but, from 
the strain amplitude h a measure of the functional combination m 5 / 3 /d, and from the polarization a measure of the 
orbital inclination i. In principle, the orientation of the gravitational wave polarization axes on the sky also informs 
us of the position angle of the orbital axis on the plane of the sky. Models for the formation of close double white 
dwarfs suggest that chirp masses m span only a relatively modest dynamical range: Webbink and Han [10] quote 
(m) = 0.42 Mq, with standard deviation a,„ = 0.14 M Q for model 4 in the population synthesis study by Han [3]. The 
source amplitude h, combined with its polarization h x /h + enables one to make a very crude estimate of the source 
distance d. Unfortunately, the gravitational wave signal alone does not permit extraction of the individual binary 
component masses Mi andM2. 

Consider a fairly typical close double white dwarf of orbital period P or b = 200 s and chirp mass m = 0.42 Mq at 
a distance d = 8 kpc. Averaging over random orbital inclinations, such a system will produce a dimensionless strain 
amplitude h = 1.0 x 10~ 22 at Earth, corresponding to a signal-to-noise ratio SNR w 70 at LISA'S expected sensitivity 
at v = 10 mHz [7, 8]. If subject to orbital energy losses due only to gravitational radiation, this source evolves toward 
higher frequency at a rate 

v = — f-^-J v 11/3 =6.3x 10~ 15 s- 2 . (4) 



This rate of period change suffices to produce a phase drift of some 20 radians in one year, rapid enough to be readily 
detectable. Measurement of V then permits determination of the chirp mass m independent of the source distance. 
The source distance itself can then be determined from the gravitational wave amplitude (see Eq. 1 above), given the 
orbital inclination from polarization (Eq. 3). Indeed, it should be possible to measure v for the overwhelming majority 
of resolvable sources. 



THE APPROACH TO MASS TRANSFER 

Suppose that our typical binary consisted of equal-mass white dwarfs, Mi =M2 = 2 1 / 5 m = 0.48 M©. If the components 
of this binary were not subject to internal energy dissipation of any sort, but remained fully degenerate up to the onset 
of mass transfer, with orbital evolution driven purely by gravitational radiation, then that onset should occur after 
inspiraling for t = 1.7 x 10 4 yr, when the orbital period had decayed to f rb = 75 s (v = 27 mHz). By that time the 
dimensionless strain amplitude h should have grown by a factor of 1 .9, and the drift rate of the gravitational wave 
frequency, v, grown by a factor of 34. 

For comparison, consider also a system consisting of two 0.7 M white dwarfs, again at a distance d = 8 kpc with 
orbital period P or b = 200 s and evolving solely under gravitational radiation. Such a system would be of particular 
interest as a possible type la supernova progenitor. It will produce a gravitational wave amplitude (again averaged over 
orbital inclination i) of h = 1 .9 x 10~ 22 at v = 10 mHz, corresponding to SNR w 130. The white dwarfs in this case are 
more compact than in the previous example, and their higher masses drive more rapid evolution. After inspiraling for 
t = 9.9 x 10 3 yr, this binary encounters the onset of mass transfer at orbital period P or b = 43 s. The strain amplitude at 
that time reaches /i = 5.3x 10~ 22 (SNR w 150), and not only is V = 3.2 x 10~ 12 s~ 2 large enough to produce a phase 



shift of 3.9 rad in one week, but V = 8.2 x 1CT 22 s~ 3 has grown large enough to produce a dramatic phase drift in its 
own right (A<j> = 54 rad) over the course of one year's observation. 

The approach to tidal mass transfer and merger in the two examples given above has been highly idealized. In 
reality, additional dissipation mechanisms come into play, notably tidal heating and tidal mass transfer up to the point 
of runaway merger. The issue of tidal heating on the approach to merger was raised by Webbink and Iben [11], and 
explored in considerably more detail by Iben, Tutukov and Fedorova [5]. It would seem almost certain that detached 
close double white dwarfs rotate asynchronously up to the point where mutual tidal perturbations become significant. 
If we suppose them to rotate slowly, as isolated white dwarfs generally do, then tides must impart roughly 1 % of the 
orbital energy of the binary to rotation in order to bring the white dwarfs into synchronism. This is not enough energy to 
lift degeneracy in the white dwarf interiors, but it is capable of substantially heating them [cf. 5], thereby significantly 
influencing conditions under which carbon ignition may take place. Nor is rotational energy large enough to alter 
qualitatively the course of orbital evolution; but the very short orbital evolutionary time scale prevailing immediately 
prior to merger, and the fact that tidal dissipation is strongest near the surface of the distorted star, imply that tidal 
heating can drive the dissipating star up to relatively high luminosities. 



Tidal Heating 

Iben et al. [5] have modeled tidal heating in a close binary white dwarf approaching merger, on the ad hoc assumption 
that the (spherically-averaged) local heating rate equals the local rate of increase in rotational energy per unit mass 
needed to maintain synchronous rotation. They find that the surface luminosities of dissipative white dwarfs vary with 
orbital period as 

^«P OTb 9 , (5) 
(see their Figure 6), with final luminosities at the onset of mass transfer approximately, 

L 2 ~ 200 (M/M ) 61 L . (6) 

These results suggest that white dwarf luminosities of order 2 L to 20 L (per white dwarf component) may 
be achieved immediately prior to tidal mass transfer in the first and second examples, respectively, posed above. 
Characteristic growth time scales for these white dwarfs are of order (L2/L2) ~ — |(P rb/Porb) = 430 yr and 50 yr, 
respectively. These quantitative estimates of tidal heating depend entirely on the effective viscosity in the envelopes of 
the heated white dwarfs, and so are inherently uncertain. Microscopic viscosity is negligible (as usual in astrophysical 
settings), but shear flows could generate enough turbulent dissipation, or amplify a seed magnetic field sufficiently, to 
drive synchronization on an inspiraling time scale. Some significant degree of tidal heating would seem likely in any 
case. 



The Onset of Mass Transfer 



The onset of tidal mass transfer in an inspiraling system does not lead immediately to runaway mass transfer and 
merger, even when the binary satisfies the conditions for dynamical time scale mass transfer, but is preceded by a 
phase in which an initial trickle of mass is driven by orbital contraction, up to the point at which the mass transfer rate 
becomes comparable to the mass of the donor, divided by the gravitational wave inspiraling time scale. With some 
idealizing assumptions (most notably, completely degenerate white dwarfs undergoing conservative mass transfer, 
apart from gravitational radiation losses), it is possible to solve analytically for the relation between mass transfer rate 
and time since contact, in the limit that the total mass transferred is small compared with the mass of the donor [11]. 

Let T72 = (dlnR2/dt)M 2 be the rate of expansion (or contraction) of the donor star in the absence of mass transfer 
(772 = for cold white dwarfs); K]l = (d\nRL/dt) be the rate of change of the Roche lobe radius of the donor in 
the absence of mass transfer (in the present case, rji = —2 / T gr , where T gr is the time scale characterizing angular 
momentum loss by gravitational radiation); £2 = (<?ln/?2/c>lnM2) a d be the adiabatic radius-mass exponent for the 
lobe-filling star; and £/, = (d\\\RL/ d\\\M2)t be the dependence of donor Roche lobe radius on donor mass in the 
presence of mass transfer, but absent any orbital angular momentum losses not coupled proportionally to the mass 
transfer rate. The relation between time and mass transfer rate can than be written in dimensionless form: 
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FIGURE 1. Dimensionless mass transfer rate, | A 3 1 , as a function of dimensionless time, t* , for dynamically stable (solid line) and 
dynamically unstable (dashed line) mass transfer, for a fully-degenerate donor star. Physical scaling factors for these two solutions 
are given in the text. 

where the dimensional time t and mass transfer rate are related to the dimensionless quantities t* and X through 
the relations t = xt* and M2 = — /1A 3 . The scale coefficients, % for time and jx for mass transfer rate, are 1 : 



Here, ji' e is the electron mean molecular weight of the donor white dwarf envelope, and \i e is the electron mean 
molecular weight of its interior. Expressions for the dimensionless functions F{q) (where q = M2/M1, the ratio of 
donor to accretor mass) and ^(A^/Mch) (the ratio of donor white dwarf mass to the Chandrasekhar mass) are given 
by Webbink [9]. Typically, they are of order F w 0.015 - 0.03, E w 20 - 40. 

There exist two branches to the solution for f*(A), shown in Fig. 1. If £2 > £l> mass transfer is dynamically stable. 
This situation holds only if the donor white dwarf is considerably less massive than the accretor, and never when 
M2 > 0.67 Mq, in the limit of conservative mass transfer. On this branch, T, Ll, t*, and A are all positive. The mass 
transfer rate asymptotically approaches M2 — — M as t -> 00. If £ 2 < £ L , mass transfer is dynamically unstable. This is 
the situation for both examples here. On this solution branch, T, jU, ?*, and X are all negative. The mass transfer rate 
becomes singular at t* = 2n/V3. 

For uniform composition and equal component masses, the scaling coefficients in the two examples give above 
are x — —1.4 yr (time from onset of mass transfer to merger t — —{2n/ 'v / 3)t = 5 yr), fx = —1.3 x 10~ 4 M yr _1 
for the typical close double white dwarf (assuming equal component masses), and T = —0.29 yr (t = 1.1 yr), 11 = 
—8.3 x 10~ 4 Mq yr -1 for the putative type la supernova progenitor. In the examples at hand, the binary components 
fill their Roche lobes simultaneously in both cases, so there is no difference in gravitational potential between the 
binary components. But for even slightly different component masses, these characteristic mass transfer rates would 
drive the accretion luminosities up to their respective Eddington limits. 



1 We restore here a factor inadvertently omitted from the expression for T in Webbink and Iben [11]. 
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THE POPULATION OF CLOSE DOUBLE WHITE DWARFS ON THE VERGE OF 

MERGER 



It is clear from the above discussion that, far from lurking in anonymity at low luminosity up to the instant of merger, 
close double white dwarfs can be expected to advertise their imminent demise by gradually brightening over the 
decades or centuries leading up to merger, first by tidal dissipation in the envelopes of their donor (less massive) 
components, then by a dramatic rise in accretion luminosity over the final year to decade preceding merger 2 . 

The obvious question is whether the event rates for close binary white dwarf mergers are high enough to be 
observable. If we adopt again the population synthesis model [3] cited above, it predicts a total galactic white dwarf 
merger rate of 0.030 yr _1 . Given the estimates obtained above for the characteristic tidal heating and mass transfer 
time scales of a typical close double white dwarf system, of order a dozen systems in the Galaxy should lie within their 
final e-folding tidal heating time scale before merger, with a probability of > 20% that at least one of these systems has 
begun tidal mass transfer. The same population synthesis model yields a galactic white dwarf merger rate for super- 
Chandrasekhar mass systems of 3.6 x 10~ 3 yr _1 , consistent with estimates of the galactic SNIa rate. There is about a 
20% probability that one of these systems galaxy-wide will lie within its final e-folding tidal heating time scale before 
merger, but only a <~ 0.4% probability that a super-Chandrasekhar mass system has begun tidal mass transfer and will 
merge within the next year. 



OPTICAL IDENTIFICATION OF LISA SOURCES 

There are important astrophysical issues that could be addressed by optical observations of these close binary white 
dwarfs in their final approach to merger. The magnitude and efficiency of tidal heating is of course one of them. 
Do the component white dwarfs ever approach synchronism? To what extent are the white dwarf interiors reheated 
by tidal dissipation? The calculations published by Iben et al. [5] indicate that such reheating could be substantial, 
and could be exercise major influence on conditions under which carbon ignition occurs in merging CO white dwarfs. 
Gravitational wave observations alone provide only limited insight into this question: if v can be measured, in addition 
to v and v, then its departure from the case for pure gravitational radiation, v = ^-V 2 /v, must be laid to orbital 
energy losses (such as tidal heating, mass transfer, or mass loss) other than gravitational radiation. Of great importance 
would be determination of the mass ratios of these incipiently merging close double white dwarfs. Gravitational wave 
observations alone constrain only the chirp mass m, while the orbital period at the onset of mass transfer is dictated 
essentially by the mean density (and hence by the mass) of the donor white dwarf. Absent knowledge of the binary 
mass ratio, the dispersion in mass ratios at a given chirp mass translates into a large uncertainty in the time-to-merger 
of a given system, even late in its evolution. But with a mass ratio in hand (and certainly in the case where the secular 
growth in system luminosity is measurable), it becomes possible to reduce the uncertainty in time-to-merger to a 
human time scale, possibly to decades or even years. The chance is remote, but nonzero, that it will be possible not 
only to identify the next type la supernova, but to predict the date of its explosion (assuming prompt detonation on 
merger) within a year or so! 

The identification of optical counterparts to LISA sources presents a major challenge. LISA'S angular resolution is 
poor, of order 

7 5 1 
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[cf. 2]. Sources at v = 10 mHz, say, will require search boxes of several square degrees; and despite the relatively 
high luminosities systems on the verge of merger may reach, they are likely to remain very compact and hot, with 
only modest visual luminosities. Furthermore, if these systems share the same scale height as old stars in the galactic 
disk, those at distances exceeding one or a two kpc are likely to be heavily reddened, and observable only in the near 
infrared. 

The key to identifying close double white dwarfs will be the fact that their ephemerides are known from their LISA 
discovery observations, and that, notwithstanding their astrophysically rapid inspiral time scales, their orbital periods 
still remain are stable enough in the short run to be useful, even if the binary phase drifts significantly. These systems 



2 In reality, the duration of this latter, early mass transfer, phase has been underestimated here. The outer envelope of the donor white dwarf will 
not be electron-degenerate, as assumed above, but relatively more extended, softening the onset of tidal mass transfer. 



will be anisotropic emitters, ellipsoidal variables at the very least, and possibly much more strongly modulated if 
caught during the onset of mass transfer. Their light will be modulated with the orbital period. The key to identification 
is then to accumulate images of the candidate field of view phased with the binary period. That period in the most 
interesting cases will be so short that it is impractical to read out a large-format CCD several times per binary orbit. If 
the field of view is sufficiently uncrowded, the image can be dithered with the orbital period; if it is crowded, however, 
one may need to separate images into distinct phase bins, for example by chopping the image synchronously with the 
binary phase into separate, non-overlapping images. 

CONCLUSIONS 

The purpose of this contribution is to stimulate the groundwork for fully exploiting the expected return from LISA. 
More realistic theoretical models are needed for tidal heating (including seed magnetic fields) and for the onset of mass 
transfer, to provide templates for the observational signatures of systems on the verge of merger. On the observational 
side, instrumental capability needs to be developed to accumulate rapid, phased observations of large fields of view. 
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